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Abstract 

Let K n be the set of all n x n lower triangular (0,l)-matrices with each diagonal 
element equal to 1, L n = {Y'Y T : Y £ I\ n } and let 

c n = jnin j/+)(Z) : is the smallest eigenvalue of zj . 

The Ilmonen-Haukkanen-Merikoski conjecture (the IHM conjecture) states that 
c n is equal to the smallest eigenvalue of YqYq , where 

0 if i < j, 

1 if i = j, 

1 ~ ( 7 ) +J if * > j- 

In this paper we present a proof of this conjecture. In our proof we use an 
inequality for spectral radii of nonnegative matrices. 
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1. Introduction and Statement of the Main Theorem 

Let S = {xi, a,’ 2 ,..., x n } be a set of distinct positive integers, ( Xi,Xj) denote 
the greatest common divisor of Xi and Xj and let £ be a positive real number. 
The n x n matrices ( S) = ((xi,Xj)) and ( S e ) = (( Xi,Xj) e ) are called the GCD 
matrix and the power GCD matrix on S , respectively. In 1876, Smith [22| 
proved that if S is factor closed, then det(5) = JC =1 f>{xk), where ip is Euler’s 
totient. Since then many results on these matrices have been published in the 
literature, see e.g. B SHIM- 
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One of interesting and active area in the study of GCD matrices is their 
eigenstructure. The first results on this subject were published in the papers 
25, 18j]. Since these results are consequences of some theorems on Riesz bases 
in functional analysis, the paper of Hong and Loewy 12] can be considered as 
the first paper on the study of the eigenvalues of GCD and related matrices 
due to the number theoretical aspect of the subject. Since their pioneering 
paper many results on the subject have been published in the literature, see e.g. 

HI U, 13, 15, 3 2(J 2l|,[22j]. In that paper Hong and Loewy investigated 
the asymptotic behavior of the eigenvalues of power GCD matrices by using 
some tools of number theory. Beside their results on asymptotic behavior of 
these matrices, in the same paper Hong and Loewy introduced a constant c n 
and used it to present a lower bound for the smallest eigenvalues of power GCD 
matrices. Let K n be the set of all n x n lower triangular (0, l)-matrices which 
each diagonal element equal to 1 and let L n = {YY T : Y £ K n }. They defined 
c n depending only on n as follows: 


c„ = mm 
Z£L n 


: n^\Z) is the smallest eigenvalue of zj . (1.1) 


Then they prove that 


A P((S s ))>c n 


min {Jeixi)}, 
Ki<n 


where J e is Jordan’s generalization of Euler’s totient and A n\(S e )) is the small¬ 
est eigenvalue of the power GCD matrix ( S e ), 


see .12, Theorem 4.2]. 

In 2008, in the light of their MATLAB calculations for n = 2,3,..., 7, Ilmo- 
nen, Haukkanen and Merikoski 0 presented an interesting conjecture about 
the constant c„. 

Conjecture 1.1. [The IHM conjecture, see fldll . Conjecture 7.1] Let Yq = 
(bij) £ M n (R) be defined by 



if i < j, 
if i=j, 
if i > j- 


( 1 . 2 ) 


Then c n is equal to the smallest eigenvalue ofYoY^ . 

Numerical evidence of the IHM conjecture as follows. Recently, the first, 
the second and the fourth author of the paper have investigated that the IHM 
conjecture holds for n = 7 and n = 8 with the help of a MATLAB code. 
Our MATLAB code running on a computei0 has verified the truth of the IHM 
conjecture for n = 7 in 23 minutes and for n = 8 in 3.5 days. Since |Ls| = 2 28 = 
268,435,456 and |Lg| = 2 36 = 68,719,476,736, it would take about 3 years to 
verify the IHM conjecture for n = 9 with the help of our MATLAB code. To 
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overcome this difficulty about time, we write a different code in C programming 
language. We use Newton’s identities (see jl6() to obtain the characteristic 
polynomial of a matrix Z in L n and we calculate the smallest eigenvalue of Z 
by using Newton’s method (see 24]) to shorten the running time. Indeed, our C 
code running on the same computer has verified the truth of the IHM conjecture 
in 30 minutes for n = 8 and in 7 days for n = 9. Thus, we have concluded that 
Conjecture [O] holds for n = 8 and n = 9. This investigation has been presented 
by the first author of this paper in [§]. 

After obtaining enough numerical evidence that the IHM conjecture can be 
true, we get the motivation to find out a proof of it. 


Theorem 1.1 (The Main Theorem). Assume the setup above. Then the IHM 
conjecture is true. 


The strategy of the proof as follows. We prove that for any matrix Y in K n , 
|Y _1 | < |Y q - I, where the matrix lb given by (11.21) . Secondly, we show that 
\Z~ l \ < |(lolo r ) -1 | for all Z £ L n . Then, by using an inequality for spectral 
radii of nonnegative matrices, we obtain a proof of the IHM conjecture. We 
conclude the paper with a discussion about further studies on the constant c n 
and a conjecture on the uniqueness of the matrix Ijj- 


2. Proof of the Main Theorem 

First we present a simple fact about a particular nilpotent (0,l)-matrix which 
we use in the course of our proofs. Here we give the proof of this fact though 
one can find in the literature. 

Lemma 2.1. Let Y £ K n and N := Y — I , where I is the nxn identity matrix. 
We denote by ( N k )ij the ij—entry of the positive integer k-th power of N. Then 
we have the following 

i) ( N k )ij = 0 whenever i — j < k, 

ii) Y- 1 = I — N + N 2 -+ (-l) 71 - 1 ^"- 1 . 

Proof. Since it is clear that N is a strictly lower triangular (0,l)-matrix, the 
proof of the first claim follows from the matrix multiplication. For the proof of 
the second claim, consider 

I n - (—N) n = (/ + N)(I -N + N 2 -+ (-l) 71 " 1 ^"- 1 ). 

Since N n = 0 and Y = I + N, we have 

y- 1 = I - N + N 2 -+ (-l)"" 1 ^" 1 . 


□ 

Now we investigate the inverse of any matrix Y in K n . In the following 
lemma we obtain a recurrence relation for the entries of the inverse of Y. 
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Lemma 2.2. LetY £ K n and N := Y—I. Also, let N = (riy) andY 1 = (oy). 
Then we have the following recurrence relation for aki ■ 

( 0 if k < l, 

aid = < 1 if k = l, 

{ - Eii 1 n uau if k>l. 

Proof. When we multiply both sides of the equality in Lemma 12. II (ii) from the 
left by —IV, we have 

-NY- 1 = -N + N 2 -+ (-l) n " 1 lV ri - 1 + {-l) n N n . 

Since N n = 0, one can easily obtain 

I-NY~ 1 =Y~ 1 . (2.1) 

By Lemma 12. II it is clear that aki = 0 if k < l and aki = 1 if k = l. Now, from 
EH), we have 

n 

aki 'y '] nuau 

t=i 

and hence, by Lemma 12.II we obtain 


fc-i 

aki = 'y \ vikiau 
i—l 


for all k > l. 


□ 


In the following theorem, we find the largest absolute value of each in 
terms of Fibonacci numbers which is a surprising result. 

Theorem 2.1. Let Y £ K n and Y~ x = (a, (? ). Then, for 1 < j < i < n, we 
have \aij\ < Fi-j, where Fi~j is the (i — j)-th Fibonacci number. 


Proof. Let N = () be as in Lemma E721 Also, let 1 < j < i < n and t = i— j. 
We will prove that | < F t by induction on t. Let t = 1. By Lemma [2.21 

we have \aj+±j\ = rij+i,j, where Tij+ij can be 0 or 1. Thus, \aj+±j\ < 1 = F\. 
Now assume that \aj+t,j\ < Ft for each t = 1,2,..., k — 1. By Lemma \2 . 2 1 we 
have aj+k,j = y nj+k,iaij. 


Case 1. Assume nj+kj+k- 1 = 0. Then \aj + k,j \ = 2 rij+k,i 


Also, 


by Lemma 12.21 aj+k-i,j = 2 n j+k-i,iaij. Since both of nj+k,i and 

rij+k-i,i for each * = j, j + 1, ■ • •, j + k — 2 can arbitrarily be 0 or 1, it is clear 
that aj+k,j and aj+k-i.j have the same values. Therefore, by the induction 
hypothesis, we obtain \aj+k,j\ < Fk-\ < Ffc. 

Case 2. Assume nj+kj+k -1 = 1- 
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Subcase i. Assume rij+k-i,j+k -2 = 0. By Lemma 12.21 we have 


\ a j+k,j\ < 


j+k-2 

i=j 


+ \a,j + k~\,j \ ■ 


Also, by Lemma [2.21 it is clear that aj+k-i,j = 2 n j+k-i,iO'ij■ Since 

both of rij + k,i and nj+k-l,i for each * = j, j +1, ■ • • ,j + k — 2 can arbitrarily be 0 
or 1, it is clear that 2 n j+k,i a ij an d dj+k-i,j have the same values. Thus, 


by the induction hypothesis, we obtain 


j+k-2 


n= 


j ^j+k : i&ij 
~\j+k —3 t 


< Fk~ i- Beside this, 


by our assumption in Subcase i, \aj+k-i,j\ = J2i=j n j+k-i,iaij ■ Since both 
of rij+k-i,i and rij + k- 2 .i for each i = j, j + 1,..., j + k — 3 can arbitrarily be 0 
or 1, it is obvious that Y^i=j 3 n j+k-i,iO-ij and a,j+k- 2 ,j have the same values. 
By the induction hypothesis, we obtain \a,j+k-i,j\ < Fk- 2 - Thus, \a,j+k,j\ < 


Fk-i + Fk -2 = Fk. 

Subcase ii. Assume rij+k-i,j+k -2 = 1- By Lemma 12.21 we have 



j+k -3 


j+k -1 

\ a j+k,j\ < 

} ^ n j+k,idij 

+ 

^ ^ Tlj+kii^ij 


i=j 


i=j+k —2 


Since 3 n j+k,i^ij and aj + k- 2 ,j have the same values, by the induction 


hypothesis, we have 


j+k —3 

i—j n j+k,iO-ij 




< Fk-2- 


In addition to this, since 


fij+k- 1 ,j+k -2 = 1 we obtain 


j+k —1 j+k —3 

^ ^ r fl'j+k,iQ j ij = j+k,j+k —2 '^-)® J j+k—2,j ^ ^ ^j+k—l^i^ij- 

i=j+k-2 i—j 

Here each (1 — Uj+k,j+k- 2 ), n j+k-i,j, ■ ■ ■ ,rij+k-i,j+k -3 can arbitrarily be 0 or 
1. Thus, ^i=j+fc _2 n j+k,idij and aj+k-i,j have the same values and hence, by 


the induction hypothesis, 


E j+k—1 

i-j+k -2 n j+k,iO-i 


< Fk- 1 . Therefore, we obtain 


\aj+k,j\ < F k -2. + F k -1 = F k . 

The principle of induction completes the proof. 


□ 


Let A = ( a,ij ), B = ( bij) £ M n ( R.), that is, the set of all n x n real matrices. 
We write A > 0 (or > 0) if all a t j > 0 (or > 0). Also, we write A > B (or > 0) 
if A — B > 0 (or > 0). In addition to this, we define |A| = (|a*y |), that is, |A| is 
the element-wise absolute value of A. The largest eigenvalue of A in modulus is 
denoted by p(A) and called the spectral radius of A. Now we fix the notation 
for the rest of the paper. 

Theorem 2.2. Let Yq = () be as in il.tfy and let Zq := YqYq ■ For all 
Z £ L n , we have \Z~ X \ < \Zq 1 \. 
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Proof. First we obtain the inverse of Yq. Let Nq = Yq — I and Nq = ( rriij ). 
Then 

[ 0 if i < j, 

y | — otherwise. 

We claim that the inverse of Yo is the n x n matrix (dj), where 

( 0 if i<j, 

Cij = < 1 if i = j, 

{ ( D' ■'/••, , if * > j. 

Since Yo £ K n by Lemma 12.21 it is clear that dj = 0 if i < j and dj = 1 if 
i = j. Also, by Lemma T2. 2 1 we have 


i—1 

Cij — ^ ^ rriikCkj 

k—j 

for i > j. Now we prove that dj = (—1 ) l ~ J Fi-j whenever i > j by induction 
on t = i — j . For t = 1, 


c j+i,j ~ m j+ij ~ 1 — 

Assume that Cj+t,j = (—1 )*Pt for all t = 1,2,..., k — 1. Recall that 


j+k-l 

c j+k,j = - E Wlj-\-k,sCsj • 

s=j 

By the induction hypothesis, if k is even then we have 

k/2 

Cy - k.jj — ^ ^ P2s—1 — Pk 
s=1 


and if k is odd then 


( fc — 1)/2 

Cj+k.j = 1 ^ ' P2s ~ Pk- 

s=l 


Thus, Cj +k ,j = (-1 ) k F k . 

Secondly, we calculate the inverse of Z 0 . Since Z 0 = YoF 0 T and Y^ - = (dj), 
we have 

n n 

(Zo 1 )a = Y, c li = 1 + E F k-i 

fe=1 k=i +1 
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for alii = 1,2,..., n. Now let 1 < i < j < n. Then 


( + )ij 


CtiCtj 

t —1 


‘"J* 


E 

t=j +1 


CtiCtj 


= (-i + J2 (-i r^Ft-iFt-j 


t=3 +1 


= ( —+ E 

4 = 1+1 


Since Z 0 1 is symmetric, for all 1 < j < i < n, 

n 

(Zo'h = + E 

£= 2+1 

Now we prove the claim of the theorem. For each Z E L n , there exists a matrix 
Y in K n such that Z = YY T . Let Y~ x = (a+. Then, by Lemma 12.21 and 
Theorem 12.11 we have 


\{Z 


-l \ 


k =1 


= Ei° fci i 2 

fc=i 

n 

= E>«i 2 

k=i 

n 

< 1+ E 


fc=2+l 
-1\ 


= (V)«| 


for alH = 1,2,... , n. 


7 




Let j > i. By Lemma \2 . 2 1 and Theorem 12. II we have 




E 


O'tiO'tj 


< 


Emm 

t =1 

n 

a ji\ + E MM 


< 


t=J+l 

n 

Fj-i + E Ft-iFt-j 

t=j +1 
| (-^0 )b | 


Finally, since Z 1 and Z 0 1 are symmetric |Z 1 | < |Z 0 | for all Z £ L n 


□ 


The following lemma is crucial in proof of the Main Theorem. 

Lemma 2.3 (See Q, Theorem 8.1.18). Let A,B £ M n (M). If |A| < B, then 
p{ A ) < p{\A\) < p{B). 

Finally, we are ready to give proof of the Main Theorem. 

Proof of the Main Theorem. Let Zq be as in Theorem 12.21 First we prove 
that the matrices Zq 1 and \Zq X \ have the same characteristic polynomial. By 
the definition of the trace of a square matrix, it is clear that 

n 

trace((Z 0 ) )= (Z 0 )i 1 i 2 ... (Z 0 )i k _ 1 i k (Z 0 )i k i 1 

ii,...,ifc=l 

for each k = 1,2, Also, from the formulae for (Zjf )ij in the proof of 

Theorem l2.2l one can easily show that sgn(ZQ )ij = (—1 for all 1 < i,j < n. 
Thus, we have 


Sgn((Z Q )iji 2 . . . (Zq (Zq )i k i 1 ) — 1 


and hence trace^Zff 1 ^) = trace^Zff 1 ) 1 *). By Newton’s identities 
obtain that Zq 1 and l-^ 1 ! have the same characteristic polynomial. 
p(\Zq 1 \) = p(Zq 1 ). From Theorem 12.21 and Lemma T2. 3 1 now we obtain 


16 ], 


we 


Thus, 


p(Z-')<p(\Z- 1 \)<p(\Zq 1 \) = p(Zq 1 ). 


Since all Z in L n are positive definite, the smallest eigenvalue of Zq is less than 
or equal to the smallest eigenvalue of Z for all Z in L n . □ 





3. Some Results and Open Problems on the Constant c n 

In the literature there are not so many results on estimating the value of 
c n . Recently, Mattila [19J has presented a lower bound for c n . Indeed, he 
proved that c n is bounded below by ( Ti 4 +9K /^ 2t ^ +lt ) • Then he showed that 

this lower bound can be replaced with ( ni+sen^+isn ) ~ 2 ~ when n is even, and 

( n 4 + 50 n^+ 48 n— 51 1 when 7i is odd. Recently, beside Mattila’s results, Altmi§ik 
and Buyiikkose |5J have obtained a lower bound for the smallest eigenvalue t n 
of the n x n matrix E^E n , where the ij — entry of E n is 1 if j\i and 0 otherwise, 
i.e., t n > (n J2k=i M 2 (^)) ■ Indeed, this bound can be used instead of lower 

bounds including c n for the smallest eigenvalues of GCD and related matrices 
defined on S = {1, 2,..., n} in the literature, see 0, UH EE EH HI] ■ After above 
studies on estimating the value of c n , we naturally raise the following problem. 

Problem 3.1. Can one improve the lower bounds mentioned above for c n ? 

On the other hand, in our investigation Q, we cannot find any matrix Y 
other than Yq in I\ n such that c n is equal to the smallest eigenvalue of YY T 
for each n = 2,3,...,9. After this observation we can present the following 
conjecture. 

Conjecture 3.1. Let n > 2. There is a unique matrix Y in K n such that 
c n is equal to the smallest eigenvalue ofYY T . In other words, if the smallest 
eigenvalue of YY T is equal to c n then Y = Yo, where To is defined by El- 
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